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ABSTRACT 

We investigate dynamical coupling timescales of a neutron star's superfluid core, taking 
into account the interactions of quantized neutron vortices with quantized flux lines of 
the proton superconductor in addition to the previously considered scattering of the 
charged components against the spontaneous magnetization of the neutron vortex line. 
We compare the cases where vortex motion is constrained in different ways by the array 
of magnetic fltix tubes associated with superconducting protons. This includes absolute 
pinning to and creep across a uniform array of flux lines. The effect of a toroidal 
arrangement of flux lines is also considered. The inclusion of a uniform array of flux 
tubes in the neutron star core significantly decreases the timescale of coupling between 
the neutron and proton fluid constituents in all cases. For the toroidal component, 
creep response similar to that of the inner crust superfluid is possible. 



1 INTRODUCTION 

Systems of interacting fermions are in superfluid pliases at low enougli temperatures. Tlius tlie core and inner crust regions 
of neutron stars consist mainly of superfluid neutrons. The neutron superfluid in the core, carrying the bulk of the mass and 
moment of inertia of the star, is the most important agent in the star's rotational dynamics. The observed timing properties, 
the period and spindown rate of the outer crust, are determined by the external magnetospheric torques as well as the internal 
torques reflecting angular momentum exchange with the core and inner crust. The neutron superfluid is coupled to the crust 
through the interactions of its quantized vortex lines with the charged components, the protons and the electrons, which are 
tightly coupled to each other by electromagnetic forces and to the crust through electron viscosity. 

The protons in the neutron star core are likely to be in a Type II superconducting phase with a uniform density array 
of quantized flux lines carrying the magnetic flux through the core in a state of minimum free energy. Although it makes up 
only a few percent of the neutron star's mass and moment of inertia, the proton superconductor plays a crucial role in the 
rotational dynamics of the neutron star by defining the physical properties and interactions of neutron vortex lines. A primary 
effect of the coexistence of two interacting superfluids, the neutron superfluid and the proton superconductor is that the two 
superfluid currents drag each other such that the neutron superfluid circulation around a quantized neutron vortex line drags 
a proton supercurrent circulating the vortex line. The induced "spontaneous" magnetization of the vortex lines make them 
scatter any relative electron and proton current strongly. This effect leads to very short coupling times between the superfluid 
core and the observed outer crust of a neutron star (Alpar, Langer & Sauls 1984). 

Observational clues and constraints on the dynamical properties of neutron star cores come from the study of pulsar 
timing, pulsar glitches in particular. Glitches are observed sudden increases in the rotation frequency. There is also a jump in 
the magnitude of the angular frequency derivative, the spindown rate of the pulsar. There are no signs of change in the external 
magnetospheric torque, as would be reflected as changes in the electromagnetic signature of the pulsar concurrent with the 
glitch. This in turn is assumed to indicate that the spin-up of the crust is due to an exchange of angular momentum within the 
star. The timescales found in the relaxation of the system after a glitch extend to the order of days and months, timescales too 
long for normal fluid viscosity Ekman processes. This is the argument, based on observations, for the presence of superfluid 
components in neutron stars. Fractional changes in the spindown rate at glitches and during postglitch relaxation indicate 
that the superfluid component involved contains about a few percent of the star's moment of inertia. This correlates very 
well with the inertia of the neutron superfluid in the crustal region. The dynamical properties of the inner crust superfluid, 
modeled in detail, explain the observed postglitch relaxation times (Alpar et al. 1993). The implication is that superfluids in 
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the neutron star's bulk, the core, must be already coupled rigidly to the observed outer crust on timescales shorter than the 
glitch rise times, which have been resolved down to about 30 seconds in the Vela pulsar (Dodson, Lewis & McCuUoch 2007). 
Calculations of the dynamical coupling times between a neutron suporfluid in the core and the crust and charged components 
of the star, based on the scattering of charges against spontaneously magnetized neutron vortex lines, indeed yield such short 
relaxation times (Alpar, Langer & Sauls 1984, Alpar & Sauls 1988, Andersson, Sidery & Comer 2006, Sidery 2009). 

The vortex lines in the core suporfluid should also be interacting with the flux lines carried by the proton superconductor. 
Since vortex lines are aligned with the rotation axis and flux lines along the magnetic axis of the neutron star, it is unavoidable 
that vortex line motions, in particular motions radially outward from the rotation axis that should take place if the core 
neutron suporfluid is to spin down, will encounter flux lines. Vortex lino flux lino junctions are energetically favorable, as the 
condensation energy cost of the normal matter phase in the central regions of both types of lines will be reduced by utilizing 
the overlap volume at a junction. This potential pinning of vortex lines to flux lines, and modes of vortex line motion against 
flux line pinning must clearly play a role in determining the rotational dynamics of the neutron star core, as was first pointed 
out by Sauls (1989) as a possible agent for pulsar glitches. The vortex line flux line interaction was also proposed as an agent 
for correlated spindown and magnetic field decay on evolutionary timescales for neutron stars in binaries, evolving towards 
the low mass X-ray binary and subsequent millisecond pulsar configuration (Srinivasan et al 1990). 

When investigating the internal dynamics of neutron stars the various interactions of the two superfiuid constituents may 
be of comparable importance. In this paper we investigate the dynamics of the neutron superfiuid in the star's core by taking 
into account, for the first time, the vortex linc-flxix line interaction together with the previously investigated effect of scattering 
of charged particles against the spontaneously magnetized vortex lines. We treat the neutron star core as a two component 
system, with superfiuid neutrons as one component and the charged particles, electrons and superconducting protons that 
are electromagnetically rigidly coupled together as the second component. If the charged fluid is treated in similar fashion as 
the normal component in a terrestrial superfluid, then comparisons can be made to the dynamics of vortex lines in superfluid 
Helium. In contrast to Helium the second component in the neutron star contains magnetic flux tubes. The interactions 
between the neutron vortices and the magnetic flux tubes will complicate the dynamics. 

There will be several possible scenarios. The first is that the energy needed for a neutron vortex to overcome the favourable 
conditions at a junction and pass through a fiux line will be so large that a vortex line effectively never passes through a 
fiux line. When the neutron vortices' motion is thus locked to another component of the star vortices are said to be pinned. 
Because we are considering pinning on fiux lines - prohibition of motion through fiux lines - there is still freedom to move 
in the direction parallel to the magnetic fiux lines. This is different from the situation often considered in the inner crust 
superfiuid where the neutron vortices pin to the crustal lattice. In this case the vortices are pinned to specific fixed nuclei 
and so the vortex has no freedom to move in any direction unless it unpins. At the other extreme the energy barrier of the 
obstruction could be so low that the neutron vortices simply pass through the fiux tubes with no hindrance. In between 
these two cases is that in which the vortices may gain enough energy to pass through via thermal fiuctuations, supporting 
some average rate of vortex motion in the direction perpendicular to the fiux lines. Thermally activated vortex motion in the 
presence of pinning barriers is known as vortex creep (Alpar et al 1984a). 

In Section 2 we review the dynamics in the presence of drag forces on neutron vortices alone, setting up the description 
of vortex motion in response to local forces as well as the ensuing macroscopic dynamics. In Section 3, we investigate the case 
of magnetic flux lines acting as absolute barriers to obstruct vortex motion. The case of vortex creep across the flux lines 
is treated in Section 4. In both cases, vortex motion in the direction parallel to the flux lines is governed by the damping 
force due to the scattering of charged particles against the vortices. We consider the creep of vortex lines through a toroidal 
distribution of flux lines in Section 5. 

We use a simple model in which the superfiuid constituents in the star, as well as the normal electron fluid rotate rigidly, 
and the rotation axes of all the constituents and the external torques on the star are aligned. The magnetic flux lines corotate 
with the charged component (the electrons + the superconducting protons). The present paper addresses the dynamics of 
spindown (or spinup) of the neutron star under external torques, leaving the case of precession to a subsequent paper. 



2 DYNAMICS WITH ONLY DRAG FORCES ON NEUTRON VORTEX LINES 
2.1 Local forces 

Before considering the global dynamical picture one needs to calculate the motion of a single neutron vortex under various 
forces. The equation of motion for a vortex line under a force F per unit length is the Magnus equation 

F = PnK X (V„ - Vi) . (1) 

The kinematic right hand side of this equation is sometimes called the "Magnus force" . Here p„ is the mass density of the 
superfluid neutrons, k is the vorticity vector directed along the vortex line, parallel to the rotation axis, with magnitude 
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K = h/2m„ , the quantum of vorticity in terms of the Planck constant and the bare neutron mass. v„ and vl are the 
velocities of the neutron superfluid and the vortex line respectively. 

The physical drag force on unit length of vortex line, arising from the scattering of charged particles against the magne- 
tization of the neutron vortex lines has the form: 

F = C (v, - Vi) = pcu-^r^)-^ (v, - vl) . (2) 

The drag coefficient C is expressed in terms of pc, the mass density of the charged component, n„ = 257„/k, the number of 
vortex lines per unit area in the neutron superfluid rotating at the (uniform) rate f2„, and r„ a relaxation time describing the 
vortex line-charged component interaction. The velocity of the charged component is denoted by Vc . 
Substituting this drag force in the Magnus Equation, we obtain: 

C (vc - Vl) = PnH X (v„ - Vl) . (3) 

Eliminating vl we obtain the force in terms of the velocity lag (v„ — Vc): 

F = p„Kf3'k X (v„ — Vc) + p„iif3k X kx (v„ — Vc) 

= pnHP'k X (V„ - Vc) - pnHP{Vn - Vc) (4) 

where the last equation holds when, as in the present model, (v„ — Vc) is perpendicular to k. The coefficients /3 and /3' are 
introduced to make contact with the conventional notation, employed in work on superfluid helium (ref.), to describe the 
strength of the dissipative (antiparallel to v„) and inertial (along ±k x t)„) parts of the force per unit length of a vortex line. 
In terms of C and p„K these coefficients are 



C/p„K I _ {C/pnK, 



2 



13= , (3'= ^"^^""^ (5) 

Superfluids with Cstrong/ PnH » 1, P « 1, /3' = 1 are said to be in the "strong coupling" regime, while Cweak/ Pnf^ « 
1, /3' = << 1 are in the weak coupling regime. 

In cylindrical coordinates with k and the angular velocities of the neutron superfluid as well as the proton-electron 
component along the z-axis, the superfluid velocity is = r^n, in the azimuthal direction. We denote the radial and 
azimuthal components of the vortex line velocity vl with Vr and respectively. All vortex lines at a given radius will have 
the same velocity vl, which does not depend on (f). Eliminating in the radial and axial components of Eq.([3|, we obtain: 



-1 



C Pnhi 



-1 „ «2 



(f2„ - f2c)r = ^ ^ — {Qn - n^y. (6) 



pc Pn ^n^v 

In this paper we shall use the results of Alpar, Langer & Sauls (1984) who calculated the microscopic coupling time r„ 
for the scattering of electrons (and the protons which are tightly coupled to the electrons by electromagnetic interactions) by 
neutron vortex lines. The scattering is due to the very strong localized magnetic field structure carried by the neutron vortex 
lines. This "spontaneous magnetization" arises from proton supercurrents circulating the vortex line, dragged along with the 
usual neutron supercurrent circulation. This mechanism, a two-superfluid effect, gives an extremely short coupling time 

.„.6.2P(pc,.)-/«(g.)^(l-^)^^^(10 - 200)P (7) 

where P is the neutron star's rotation period, Pc,i4 is Pc in units of 10^* gm cm~^, rUp and Snip are respectively the proton 
mass in vacuum and its deviation from the vacuum value in the neutron star medium, with the effective proton mass given 
by rrip* = rUp — Snip (Alpar, Langer & Sauls 1984, Alpar & Sauls 1988). The prefactor (10-200) describes the variation of 
r„ with density, proton fraction and proton effective mass throughout the neutron star core, with shorter coupling times at 
higher densities. The proton effective mass is not known very accurately, theoretical calculations giving rrip* /rUp ~ 0.5 — 0.9. 
We have used the results of (Chamel & Haensel 2006, Baldo & Ducoin 2009) at nuclear matter density pnuc — 2.8 x 10^'' gm 
cm"'' and at 3p„uc. The range (10-200) P also reflects the variations in the calculated parameters as reported in these two 
papers. Eq. Q gives 

since the charged component makes up only a small fraction of the density in the neutron star core, pc/pn ~ 5 x 10~^. 
Thus the charge scattering by spontaneous magnetization of the neutron vortex lines leads to short coupling times, which are 
nevertheless in the weak coupling regime in the sense of Eq. ||8|. The coefficient in Eq. Q simplifies to 

Vr = — (!^„ - f2c)r = — (!^„ - Qc)r. (9) 

Pn ^^^n'^v Pnl^ 

The neutron star core superfiuids would be in the strong coupling regime if much stronger coupling mechanisms existed 
between the charges and the neutron vortex lines (Sedrakian & Sedrakian 1995). Strong coupling would keep the neutron 
vortex lines almost in corotation with the charged component and flux lines. The presence of flux lines therefore should not 
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alter the dynamical coupling times significantly. The dynamical coupling times obtained from the particular range of drag 
coefficients Cstrong — T] oi Sedrakian & Sedrakian (1995) are much longer than the observed glitch rise times. We will derive 
general expressions for dynamical coupling times in terms of the parameters (3 and /3' which can be used for both weak and 
strong coupling. The model of Alpar, Langer & Sauls (1984), which is in the weak coupling regime, will be used for the 
applications throughout the paper. 



2.2 Macroscopic Dynamics 



We are now ready to introduce equations for fl„ and Qc, the two macrosopic dynamical variables of the rotational dynamics. 
The superfluid is able to spin up or spin down only by radial motion of the vortex lines, 

r r 

This equation, appropriate for axial symmetry about the rotation axis, follows from the fact that a superfluid's vorticity 
resides in quantized vortex lines, whose number density is conserved: 



When integrated over a cylindrical volume vrr times a unit length. Gauss' Theorem leads to 



Vrdl 



(11) 



(12) 



where the line integral is on a circle of radius r around the equator of the star in the x-y plane with the z-axis parallel to the 
axis of rotation. Eq. ( 10 1 follows with dl = rd(j> and axial symmetry. 

The last equality holds in the case of uniform rotation. Substituting the model result in Eq. ^ for Vr gives 

1 



where 





Pn P' 



(f^Ti ^c) 



~c~ 



+ 



c 



P_ 

4:71 



(13) 



(14) 



and uj = ^ln — fie denotes the lag between the rotation rates of the two components. This is complemented by the angular 
momentum evolution of a two-component star under an external torque Next, 

InCln + iMc = Next = iClac (15) 

where / = /„ + 7c is the total moment of inertia and Qoo = Next /I denotes the steady state spin down (or spin up) rate 
reached asymptotically - hence the suffix oo - when both components are spinning down (or up) at the same rate. In the 
present linear drag force model, Eq.([2|, initial conditions away from the steady state actually relax exponentially to the steady 
state. Eq.(15l can be written, at each radius r in an axisymmetric model, as 

p„Q„ + pc^c = piloo (16) 

where p = + Pc is the total density. In the neutron star core matter is mostly neutrons, with Pc/ p of the order of a few per 



cent. From Eqns. ( 13 1 



( 16 1 we obtain the solution for the lag in rotation rate between the two components 
uj{t) 



where the dynamical relaxation time r is 



and the steady state lag is given by 



Pc 

P 



pnK 
C 



+ 



c 

Put 



LOoo] exp(— t/r) + ujo 



P_ 



P 



I i i oo I 'v 

pc 



(17) 



(18) 



(19) 



The solutions for f2„(t) and Q,c(t) for any initial conditions can be obtained from the solution for the lag, Eq. (171, and 
the integration of Eq. (161. These solutions are also characterized by exponential relaxation with the dynamical relaxation 
time T. We note that the result of Alpar & Sauls (1988) for the dynamical relaxation time, r = p„ / pcT^, was an artefact of 
keeping Q,c constant (Andersson, Sidery & Comer 2006). 

The physics of any particular model of superfluid - normal matter coupling is described by an equation like Eq.( 12 1. When 
the dependence on the lag lo is linear the coupling is described by r„. The dynamical relaxation time r and the steady state 
lag LOao scale with r„, as given in Eqs. (18 1 and (19 1. These generic relations hold in two-component models with different 
physics, as long as the coupling is linear like in Eq.(13l and angular momentum evolution is given by Eq.(15l. 
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Figure 1. A sketch of the neutron star core. In the frame fixed with the charged fluid the magnetic vortices are stationary and parallel 
to the a-axis. The neutron vortices can move in the x direction, encountering only the drag force. 



3 DYNAMICS WITH VORTEX LINE - FLUX LINE PINNING AND DRAG FORCES 
3.1 Local forces 

Let us now introduce flux lines. The proton flux lines act as barriers offering obstruction to neutron vortex motion in the 
direction perpendicular to the flux lines. To impose this we introduce a cartesian coordinate system in which the z-axis is 



parallel to the axis of rotation (See Figure (3.1 1). We fix the a;-axis in the direction of the proton vortices (projected on the 



x-y plane) and so the coordinates will rotate with the proton fiuid. The neutron vortices can move in the x direction, moving 
parallel to the flux lines, but are obstructed in the y direction, perpendicular to the flux lines. 

To include the obstruction by the flux lines we use a directionally dependent drag coefficient and rewrite the balance of 
forces as 



Cx (Vx^c - Vx,l) X + Cy {Vy,c - Vy^L) f = Pn K X (V„ - Vl) 



(20) 



Noting that this stiU holds in a rotating frame we can set Vx,c = Vy^c = by using the coordinate system fixed to the proton 
flux hnes, proton superconductor and electrons (the charged, c, component). The x and y components of this equation are 



Cy'^y,L — pnl^ {f^x,n 1^x,l) ■ 



(21) 
(22) 



The case of absolute pinning in the y direction corresponds to the limit Cy ^ oo . Equation ( 22 I shows that for a reasonable 
solution we must have Vy^L = 0. To put it more carefully, the physical force in the y direction is not necessarily a linear drag 



force, it is some probably highly nonlinear force that imposes Vy^L = 0, so that Eq. (22 1 becomes 



Fy 



(23) 



As the vortices are not constrained by the fiux lines in the x direction Cx is simply the electron scattering drag coefficient C. 

(24) 



Substituting this into Eq.(21 1 gives us the force on unit length of vortex line as 

Fx 

In the y direction we have, 



-Cv, 



x,L 



pn f^Vy.n 



Fy = P„liVx 



^ f^y .n — ■ pn ^"^x^n Pn ^y^n ■ 



(25) 



In the inertial frame where Vc 7^ the force per unit length is: 

F„ = p„K X (v„ - V, 

= p„K X (v„ - V, 



2 2 
PnK 



C 



(,v„ • y)y 



P^i^-pji^ri ■ y)y 



(26) 



which clearly does not have the symmetries of Eq.(|4|. 

Assuming an array of vortices of constant density n„ and expressing the velocities in terms of the angular velocity lag lj, 
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we find the force density on the neutron superfluid to be 

Fx,n = —Pninriv) LUX 



P — 



(27) 
(28) 



We should also note that the form of the force in ( 24 1 and (251 in the limit of strong damping by electron scattering 
(C — > oo) is the same as in the case with no magnetic flux lines. In this strong coupling limit the vortices are locked into the 
charged component (proton superconductor plus normal electrons) whether the magnetic flux lines are present or not and so 
the vortex dynamics are not modified by the fiux lines. 



3.2 Macroscopic Dynamics 



Up to this point we have derived the macroscopic dynamics from considerations of vortex currents, as in Eq.( |12[ ). The dynamics 
of our two component system can also be formulated in terms of torques. In this section, for illustrative purposes, we shall 
present both derivations. 

The equations of motion in a frame rotating with the proton fluid will be given by 



Ix^x + Ix^c X 17j 



r X VxdV 



star 



where X denotes either of n and c. 

For the neutron fluid component we find 



fsta 



I„fl„ + I„flc X f2„ = —2p„Q„uj I ^xrxx + rxy\^y + —x J J- dV 



Most of the integral vanishes by symmetry so we are left with 
For a spherical star of constant density 



jjx'zdV 

tar H 



.„ / x^dV = \pr.VR^ - ^" 

J Star 



(29) 



(30) 



(31) 



(32) 



5" 2 

where In is the total moment of inertia of the neutron fluid. Taking the angular velocity vectors initially in the z direction, 
only the magnitudes of fix will change, according to 

f3 



(33) 
(34) 



Alternatively, Equation (33 1 can be derived using Eq.{ 12 I, the relation between spindown rate and vortex current. In the 
present case 

(35) 



Vr = V:c,L COs((j!>) = -^Vy,n COs((/!>) = -^VLO COS {(j)) , 



where we have used Eq.(24l, and substituted Vy,„ = rcJC0s((;6). The integration of radially outward vortex current leads to 
Eq. (33 1, in the form 

with 



C 



-Q„uj 



C 



In the presence of an external torque Next aligned with the angular velocity vectors, Eq. ( 34 1 becomes 

Ic^c — ~~ In^n ~\~ Next-, 

SO that the equation for the lag is 

^ P I Next 



(36) 
(37) 
(38) 
(39) 



f3' Ic Ic ■ 

Observed rotation rates of neutron stars have variations of at most a fraction AQ,c/^c ^ 10~^ even in the largest glitches. 
The lag u must therefore be small compared to Qc and Q,„, and (39 1 can be solved with f2„ set to a constant representative 
value rio- The solution is 



io{t) = (t<j(0) — LOao) exp ( ^ j 



(40) 
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which is just Eq.(17l with the relaxation time rpi„ appropriate for the present situation. The dynamical timescale rpi„ over 

(41) 



which the two fluids couple is related to the timescale r„,pi„, with the generic relation 



leading to 



Pc 



C 



(42) 



and Woo ~ |f^oo|TVi,pin- The solutions for Q„{t) and Qc{t) are obtained from Eqs. (38l and (401. Thus the introduction of flux 
line - vortex line pinning as an absolute barrier to vortex motion in one direction reduces the coupling time by a factor 

2 



^2(3' ^2 



(—] -0(3.2x10^^-8x10"'°). 



(43) 



4 VORTEX CREEP 

4.1 Vortex Creep in the Axially Symmetric Case 

At finite temperature vortex lines can pass (creep) through flux tubes via thermal fluctuations. We shall model the resulting 
vortex creep velocity in the y direction by generalizing the model developed for the neutron star crust superfluid where vortices 
creep against an isotropic distribution of point pinning centres (Alpar et al 1984a). 

We flrst review this familiar model to set the stage for anisotropic creep against flux lines. For a vortex line to unpin 
there is an energy barrier of size Ep to overcome. The rate of spindown of a superfluid is determined by the vortex velocity 
in the radially outward direction from the rotation axis, as we saw in Eq. ( |10[ ). The creep velocity in the radial direction is 

Vr = vo exp I -y^ — ) — exp I TTFr~ I (44) 



kT J ' \ kT 



where vq is the average trial velocity against pinning potentials, Ep-out and Ep-in are the pinning energy barriers for radially 
outward and inward motion, respectively. If there is no lag between the two fluid constituents, uj = 0, then there is no difference 
in pinning energy barriers for vortices creeping outwards or inwards; Ep-out ~ Ep-in and so Vr — 0. When the star is spinning 
down, driven by an external torque, the crust and charged constituent is rotating faster than the neutron superfluid {uj > 0). 
Vortex motion in the outwards direction will reduce this lag and so is energetically favourable. 

A pinned vortex line is dragged along with the crust and charged fluid. A velocity difference will be set up between 
the line and the neutron superfluid. According to the Magnus equation, Eq.|T|, a radially inward force must be applied on 
the vortex line by the pinning centres. The pinning energy barrier in the outward direction is reduced by an amount AE 
proportional to the lag lu, while the barrier to motions in the radially inward direction is enhanced by the same amount. The 
pinning barrier Ep can sustain a maximum lag Ucrit', there is not enough pinning force to keep vortices pinned at higher lags: 

Ep = b^Fcritical = b^pnKLUcritr (45) 

where b is the distance between pinning sites and f , the coherence length, is the radius of the vortex core, the distance scale 
over which the pinning force acts. The second equation relates the critical pinning force per unit length of vortex line to the 
critical lag through the Magnus equation. At the actual lag u < ujcrit the correction AE to the pinning energy Ep is 

AE = b^F = b^pnKur. (46) 
The energy barrier in the outwards direction is given by 

Ep-out = Ep-AE (47) 

= Ep(l- . (48) 

Vortex motion in the radially inwards direction faces an energy barrier enhanced by A_E, so that 

Ep-,n = Ep(\^ —\ . (49) 

\ '^crit I 



This finally gives the radial creep velocity 



2.oexp(-§ sinh (§£-). (50) 



Substituting the radial creep velocity in the equation for superfiuid spindown, Eq. (lOl completes the creep model for the case 



of an isotropic distribution of point pinning centres. There will be no creep in the azimuthal direction. 
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4.2 Vortex Creep Across Flux Lines 



In the presence of flux lines, vortex motion in the y direction will be determined by creep through the flux lines. We model 
vortex motion in the y direction in terms of creep alone, VL,y = «creep('^), ignoring the effects of drag forces. This approach 
is appropriate in the case of weak coupling drag forces, and is justified for the drag forces given in Eqs.([7| and Q. Strictly 
speaking, a vortex line moving in the y direction experiences drag forces between encounters with flux lines. For weak enough 
drag, the governing rate (bottleneck) for motion in the y direction is given by the very slow rate of creep across the flux lines, 
so that the effect of weak coupling drag forces can be neglected. Strong coupling, with drag and creep having comparable 
effects on the motion of vortex lines in the y direction requires a more complicated treatment which will not be adressed in 
this paper. 

To flnd the appropriate expression for Vcreep we must recalculate Ep-out and Ep-i„ or rather in this case the effective 
pinning energies of the positive and negative directions along y {Ep+ and Ep- respectively). Again any bias must be due to 
vortex motion in the radially outward or inward direction. This is taken into account by introducing a factor r • y = sm{(p): 



E., 



Er, 



CO 



The vortex lines' velocity in the y direction, due to vortex creep is thus 

E^. 



sinhfl^^^^^ 



(51) 



(52) 



In the X direction, vortices move parallel to the flux lines, at velocity vl.x encountering only the drag force due to scattering 
of charged particles off the magnetized vortex core. The components of the Magnus equation Eq.(l) with anisotropy in the 



rotating frame of the charged component c, analogous to Eqs.(21 1 and (231, are 



Cvl,x — PnH {VrL,y — Vcreep) = pn^ {iOT COs((j!)) — Vcreep) (53) 

Fy{u) = p„K {vn,x ~ vl,x) = PnK {—ior sin{(j)) — Vl,x) ■ (54) 

Since the vortex line velocity in the y direction can, in principle, be a nonlinear function of u, the physical force per unit 
length Fy^creep representing the interaction of the creeping vortex with the flux lines can also be a nonlinear function of ui. 

With these two different vortex line velocities in the x and y directions, the equation of superfluid spindown, Eq. ( 12 1 
becomes 



Substitution from Eqns. ( |52[ ) and ( |53[ l leads to 

jf [Pn^ / 2 



[vl,x cos(<^) + VL,y sin(0)] ro 



j> {Ljrcos^{4>) - iicreep('i'sin(<;/>)) cos((?!))) +«creep('jjsin((/)))sin(<;/>)j rd(t>. 



(55) 



(56) 



The second term clearly integrates to zero. Evaluating the integral in the flrst term, our equation of motion is now given by 



'~cr 



2n„«o 



- exp 



kT 



sinh 



£p 



■ sm{(f)) 



sin{(l))d(f) 



(57) 



Both of these terms have the same sign, and so taking the flrst term alone will give us an upper limit on the lag ui. Taking 



the Vela pulsar as an example, we have f2 



-10" 



rad s"^, n « 70 rad s"\ The range of values of C/(p„k) quoted in Eq.Jsk 



imply that the lag will be less than ~ 6 x 10 ^® rad s ^ . This bound implies that the creep term is also linear in the 



since the argument of sinh is very small, as can be seen using Eq. (451, with a temperature T — 10 K and typical parameter 



values for the neutron star core. Eq.(57l becomes 



exp 



Ej, \ Er, 



(58) 



C R \ kT J kT LUcrit 

where we have substituted a constant flducial value Q,o for the slowly varying and taken R = lO'' cm to represent the 
radius r. Introducing a relaxation time 



R 

2Q.„vo V kT 



kTuJc 



Er, 



and an effective relaxation time 



Eq. (1581) can be written as 



(59) 
(60) 
(61) 
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As in the earlier sections, treating this coupUng together with the angular momentum evolution, Eq.(38l, leads to solutions 
of the form given in Eq. ( 40 1 with the dynamical relaxation time 

T^If = — 'Tn.e// = (62) 

where the dynamical relaxation time for creep (in the y direction) is 

"Tcreep — '^n. creep- (^^) 

p 

The effective relaxation time is the shorter of rpi„ and Tcreep- The ratio of these relaxation times is 



R f Ep\ kTuJcrit PnH^n .„.s 

exp (64) 



Tpin 2Q„vo \kT J Ep C 

We take T = lQ^K,kT « 10"^ MeV, a typical density p = 3 x 10^" g cm~^, vo = lO'^ cm s"\ R = lO'' cm, the coherence 
length 5 = 10~^^ cm and the distance between pinning junctions (flux line spacing) fe = 5 x 10~^°(-Bi2)~^''^, where B-12 denotes 
the magnetic fleld in the neutron star core in units of 10^^ G. We find that Tcreep is the longer timescale unless Ep/kT < 9. 
With estimates of the pinning energy _Bp ~ 0.1 — 1 MeV (Ruderman, Zhu & Chen 1998, Srinivasan et al 1990) and kT « 10~^ 
MeV, the creep rate is likely to be negligible, and the core superfluid remains in the very tight coupling with the charged 
component imposed by the presence of the flux lines. This result depends only logarithmically on the parameter values we 
have used. 

If the pinning energy is closer to the thermal energy {Ep/kT < 5) then the creep process determines the relaxation of 
the core, r^// = Tcreep, somewhat alleviating the almost rigid coupling to the charged component inferred in Eq.(|43[). 



5 TOROIDAL FIELD 

Recent stability simulations of the type of magnetic fields expected in neutron stars (Braithwaite 2008) suggest the presence 
of a large toroidal contribution to the internal magnetic field. In this case the fiux lines prevent radially outwards motion of 
the superfluid vortices in all directions, so that any spin down of the superfluid must be due to creep. Because of azimuthal 
symmetry, the macroscopic dynamics will in fact be the same as in the creep model for the crust superfluid. If toroidal flux 
lines coexist with the uniform array of flux lines aligned in the x direction, reconnection will lead to a toroidal field; in any case 
there is no "easy" direction with vortices moving outwards against drag forces only. We therefore assume that the toroidal 
field is confined to an equatorial belt surrounding but not overlapping with the core region with poloidal field lines (here 
poloidal means a uniform parallel array of fiux lines). The dominant physical process constraining the radial vortex current 
is creep against the torroidal flux lines. The radially outward vortex current will have its dynamics characterized by the local 
values of the lag uj. In the region of toroidal fiux lines the lag lj may have quite different values from the rest of the core. 
The equation of motion is 

^^oVo_f Ep\^_^^^^fEp 



Vortex creep can operate in the linear or nonlinear regimes (Alpar, Cheng & Pines 89), depending on a comparison of the 
spin down rate determined by the external torque and the efficiency of thermal creep, as refiected in the magnitude of the 
sinh function 

sinh ( = exp (^) = 8r.p.n.o^.^^o / Ep\ 

\kTLOcru) |n„|7? \kT) R ^\kT) ^ ' 

where Tspindomn = no/2|f2| is the characteristic spindown age of the pulsar. For the Vela pulsar parameters, the estimate of 
the sinh is less than 1 if Ep/kT < 31, or Ep < 0.3 MeV, if T = 10** °K. In this case Eq.(|65| is linearized to give 



^ Aflovo f Ep\ Ep Lo 2uj UJ 

^" " n~^^^\Tk]lJT = = • ^^^> 

^ \ '^■-^ / ^crit Tn, creep Tn^tor 

Thus, in the linear regime creep across toroidal flux lines the relaxation times T„^tor and Ttor are half the corresponding 
relaxation times r„, creep and Tcreep for creep across flux lines aligned in the x direction. This is easy to understand since creep 
operates in the r direction, that is, both x and y motions are allowed, while in the latter case creep takes place only in the y 
direction. The relaxation time of the toroidal region has a sensitive dependence on Ep/ (kT) : 

Pc Pc R f Ep \ kTucrit in-Srr f Ep\ .„0N 

nor = jT„,or = ^ ^Xp j -10 Ts CXp - j S (68) 

Thus Ttor ~ 3 X lO^s for Ep/kT = 31, the limiting value for linear creep. If Ep/kT is somewhat smaller, the dynamical 
relaxation time is reduced exponentially; for example Ttor ~< 2 x lO"** s for Ep/kT = 10. 
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6 CONCLUSIONS 

If flux lines of the proton superconductor are not taken into account, the dynamical relaxation time for the core superfluid is 

T = — r„ = — r„ ?^ r„ (10 - 200)P. (69) 
P P 

as given in Eqs. ([7| and (181, employing the results of (Alpar, Langer & Sauls 1984). The range of uncertainty corresponds 
to the uncertainty in the representative values of the effective proton mass in the neutron star core. 

We have found that the inclusion of the proton flux lines strengthens the coupling between the neutron superfluid in the 
core and the charged component of the star, leading to substantially shorter dynamical coupling times. When flux lines are 
included as absolute barriers in a poloidal like conflguration, the neutron superfluid spins down by the outward motion of the 
vortex lines in the "easy" direction, parallel to the flux lines, checked only by drag forces. The dynamical coupling time is 
reduced by the factor 2/3' ~ 2[C/ (pnn)]^, where C is the coefficient of the drag force and pnK is the Magnus inertial coefficient, 
yielding 

^-3.2x10^^-8x10-'°. (70) 

T 

Absolute pinning to poloidal flux lines forces the neutron superfluid into effective corotation with the charged component 
and outer crust of the neutron star. Allowing vortex creep across the flux lines does not effect this conclusion. For relevant 
neutron star core temperatures and pinning energy estimates for flux line - vortex line junctions, the timescale for creep across 
the flux lines is typically much longer than the timescale Tpi„ of dynamical coupling achieved by vortex motion parallel to 
the flux lines alone, so that the effective dynamical relaxation time with creep present is not changed, Te// = Tpin. Thus the 
presence of flux lines reduces dynamical coupling times drastically and clamps the neutron star core superfluid to the crust 
and charged component. 

The likely toroidal arrangement of flux lines circumscribing the core superfluid does not allow any anisotropy in vortex flow. 
All vortices have to creep across the toroidal flux lines if they are to move out radially. If Ep/kT is small enough {Ep/kT < 31 
for Vela pulsar parameters, with a similar critical value for other pulsars) creep is linear, the dynamical relaxation time 
is simply the creep relaxation time Ttor ~ 10~*r8 exp (^^^ s. The toroidal region may well contain values of the pinning 
energy with Ep/kT large enough for the nonlinear creep regime. In the nonlinear regime, corresponding, for the Vela pulsar 
parameters, to Ep/kT > 31, or Ep > 0.3 MeV, if T = lO^K, The sinh function is approximated with sinh(2;) — l/2exp(2;). 
The resulting highly nonlinear dynamical model has been successfully applied to the dynamics of the inner crust superfluid, 
and invoked to explain basic features of postglitch and interglitch response in pulsars (Alpar et al 1984b). With its location 
at the boundary between the neutron star core and superfluid inner crust, the region with toroidal flux lines could have a 
moment of inertia fraction Ip/I ~ O(10~^), comparable to the fractional moment of inertia of the superfluid inner crust. The 
steady state lag uioo — LOcrU and the relaxation times for nonlinear creep in the toroidal region may also be similar to the lag 
and relaxation times of the inner crust superfluid. The toroidal region will then have postglitch response signatures similar 
to those of the inner crust superfluid, and will be responsible for some of the nonlinear response previously attributed to the 
inner crust superfluid. This makes it impossible to draw unambiguous lower limits on the moment of inertia fraction of the 
inner crust superfluid to use for constraints on the neutron star structure and equation of state. 

Nonlinear creep implies a much larger value of the steady state lag uJao in the toroidal region than the typical lags in 
the linear regime, which we found to be extremely small in the neutron star core with its poloidal (uniform) array of flux 
lines. This change in uj means a sharp increase in f2„, and therefore a large density of neutron vortex lines in the boundary 
between the core superfluid and the toroidal region. This may well be an axially symmetric vortex reservoir to provide the 
source of glitches when vortex lines are unpinned catastrophically. Thus the toroidal region could become an alternative, or 
a complement to the inner crust superfluid as a source of glitches as well as being a prominent agent of postglitch relaxation. 
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